Imaging spin properties using spatially varying magnetic fields 
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We present a novel method to image spin properties of spintronic systems using the spatially 
confined field of a scanned micromagnetic probe, in conjunction with existing electrical or optical 
global spin detection schemes. It is thus applicable to all materials systems susceptible to either of 
those approaches. The proposed technique relies on numerical solutions to the spin diffusion equation 
in the presence of spatially varying fields to obtain the local spin response to the micromagnetic probe 
field. These solutions also provide insight into the effects of inhomogeneities on Hanle measurements. 

PACS numbers: 75.76.+j, 75.40. Gb, 85.75.-d, 75.78.-n 



Spin electronics is a rich and rapidly growing field that 
integrates magnetic and nonmagnetic materials to enable 
a variety of spin polarized transport phenomena [TJ |2] . 
Much progress has been made in understanding injection 
of spins into non-magnetic semiconductors, the manipu- 
lation of those spins within the semiconductor, and both 
local and global detection of the spin states [3H6]. Sev- 
eral recent results have brought us closer to electrically- 
controlled room-temperature spintronic devices [7HT5]. 
However, much is still not well understood, particularly 
regarding the mechanisms that degrade spin polarization, 
in the often complex magnetic environment that prevails 
in spintronic devices. Spin lifetimes inferred from mea- 
surements are often much shorter than theoretical expec- 
tations and the influence of impurities, inhomogeneities 
and other departures from ideal structures remain poorly 
understood. 

A technique which could spatially resolve the local vari- 
ations of spin density and other key parameters, such as 
spin lifetime, in multicomponent spintronic devices would 
have enormous impact. Here we show that local infor- 
mation about the spin properties can be encoded in the 
variation of the global spin polarization in response to 
a spatially scanned magnetic field from a micromagnetic 
probe. In analogy with electrostatic tip-induced scatter- 
ing in a two-dimensional electron gas [16] the local mag- 
netization is modified to an extent determined by both 
the vector field of the micromagnetic probe and the local 
spin properties. This local perturbation is detected by 
established global spin polarization measurement tech- 
niques (such as spin photoluminescence or non-local volt- 
age) which, in themselves, need not provide spatially re- 
solved information, hence the method is applicable to any 
materials system (whether optically active or inactive). 
A schematic of our proposed imaging technique is shown 
in Fig. [TJa). As a magnetic dipole is scanned above the 
sample the dependence of a globally detected spin signal 
on position of the probe provides a map of the sample's 
spin properties. 



Spin dynamics in nonmagnetic materials are given by 

<9S S 

-=D s V 2 S + 7 BxS-- + G, (1) 
at r s 

where S is the spin density, D s is the spin diffusion con- 
stant, 7 = gfiB/fo is the gyromagnetic ratio, B is the 
total magnetic field experienced by the spins, r s is the 
spin relaxation time, and G represents the spin genera- 
tion term, e.g. irradiation with circularly polarized light. 
S is a function of time t and spatial position r = (x, y) 
for the 2D sample that we consider. 

The spatially dependent response of the spin polariza- 
tion to the magnetic field of the probe tip can be calcu- 
lated by means of a numerical technique we have devel- 
oped for solving the spin diffusion equation in the pres- 
ence of spatially varying vector magnetic fields. Vector 
precession is an added complexity which is not present in 
the case of (scalar) charge deflection by a scanned gate. 
Our technique gives steady state spatial maps of spin 
density as a function of the location of the micromag- 
netic probe over the sample. The numerical solution 
is obtained from the time-dependent form of the diffu- 
sion equation using Euler's iterative method (please see 
supplementary material for further details). Though we 
do not consider electric fields nor take B, G and r s to 
be time-dependent these situations are straightforwardly 
handled by our approach. 

We consider spins injected uniformly into the sample 
by some means (uniformity is not a necessary condition). 
The injected spins and the scanned dipole are assumed 
to be oriented along z. This sample has a small localized 
region in which the spin lifetime is five times longer than 
the rest of the sample. Panels (b) and (c) of Fig. [ljshow 
spatial maps of S z , the component of the spin density 
parallel to the injected spin direction, for two positions 
of the probe relative to the lifetime inhomogeneity. Panel 
(d) presents the image of S z , g , the integral of S z over the 
entire sample, as a function of the position of dipole over 
the sample. The region of slower relaxation rate is clearly 
distinguished from the rest of the sample. 
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FIG. 1: [Color online] Calculated vector magnetization 
demonstrating imaging of local properties using global mea- 
surements: (a) Measurement geometry of our proposed ex- 
periment. Embedded in a uniform plane is a small region 
(dotted line) in which r s is five times larger than throughout 
the remaining area. The 2D sample is uniformly pumped with 
spins. A probe is positioned above the sample at a distance 
z p and its x^-coordinates are scanned to obtain the image 
presented in panel (d). (b,c) Spatial map of S z (x,y) when 
the enhanced r s region is found at (2L S ,2L S ) and (6L S ,6L S ) 
relative to the tip. (d) S z ,g, the globally integrated S z over 
the entire sample as a function of probe co-ordinates (x p ,y p ). 



The sensitivity of spin polarization to the position of 
the probe arises from the interrelationship of the two 
mechanisms that alter spin polarization: spin relaxation 
and dephasing of spins in an ensemble as they precess 
about the spatially inhomogeneous magnetic field vec- 
tors. Fields parallel and perpendicular have contrast- 
ing impacts since magnetic field parallel to the injected 
spin direction will tend to maintain this orientation and 
maximize spin polarization, while transverse components 
cause precession, thus dephasing and a reduction of spin 
polarization. This competition between the various com- 
ponents of the magnetic field is well captured by a quan- 
tity 0% = j 2 Bj_r 2 /(l + 7 2 £?yT s 2 ), which governs the av- 
erage projection of the spins along the ||-axis (e.g., for 



uniform systems, the steady-state value of S\\ is propor- 
tional to (1 + 

We are interested in spatially varying fields arising 
from a micromagnetic tip. The dipole field arising from 
a magnetic moment m = m z z located at (0, 0, z p ) over 
a sample in the x?/-plane results in a spatial variation of 
spin density S z shown in Fig. [2] (left axis) along with 
Ob (right axis). Here we assume a sample with uni- 
form r s and uniform injection of spins oriented parallel 
to z and approximate the probe field by Bdi p (x,2/) = 
^[(3R(m • R) - mR 2 ]/R 5 , where R = ix + yy + z p z 
and /in is the permeability of free space. Fig. [2] shows S z 
both with and without diffusion. The anticorrelation be- 
tween S z and Ob is particularly evident when diffusion is 
neglected because a given spin experiences a single mag- 
netic field throughout its lifetime. Diffusion smears out 
the sharp features, and if the minima in S\\ (for D s = 0) 
occur closer than a spin-diffusion length L s away from 
the central peak, then the peak will collapse leaving a 
localized volume of suppressed spin density. 

This local suppression of the spin polarization by the 
probe tip field provides the basis for scanned probe imag- 
ing of a spintronic sample: the impact of this tip field is 
determined by local sample properties. The image con- 
trast in Fig. [IJd) results from the fact that long-lived 
spins are more strongly dephased because they experi- 
ence the tip field for a longer time before being replaced 
by fresh, well-oriented spins. Thus, the globally aver- 
aged spin density of the sample will be more strongly 
suppressed when the dipole is over the region of longer 
r s . The size of this effect depends on the "global" size 
over which the spin polarization is integrated. 

To estimate the spatial resolution of our imaging tech- 
nique, we consider the effect of a linear magnetic field 
gradient: B = a(yy — zz) (see Fig. |3|. The injected 
spins are uniformly oriented along z. The resulting vari- 
ation of S\\ along the ?/-axis, for a = (7T S L S ) _1 , where 
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FIG. 2: [Color online] (a) Steady state spin density S z (x,y) 
for y = in the presence of the field of a dipole for D s — 
(red) and D s r s /z p = 1 (blue). The dipole m z z [here \i{)Vn z — 
107TZp)] is located at r = (0, 0, z p ) and uniform injection is 
assumed. The green curve shows the spin dephasing factor 
B {x)\ note that 6 B = 1 S z = 0.5 for D s = 0. 
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FIG. 3: [Color online] The resolution length scale L r as a 
function of magnetic field gradient a, where length is scaled 
by the diffusion length L s , and field is in the units of (7T S ) _1 . 
(inset) Variation of S z along the y-axis in the presence of a 
field of the form B = a(yy — zz). A Gaussian fit to the 
numerically simulated data S z is shown. The width of the 
Gaussian may be considered to be the resolution for a given 
gradient. 



(t^s) -1 is the half- width of the conventional Lorentzian 
Hanle curve [17 , is presented in the inset of Fig. [3j A 
Hanle curve is simply a plot of the dependence of spin 
polarization on a transverse magnetic field. In the ab- 
sence of diffusion, spins at different positions along the 
axis will sample different fields, resulting in a map- 
ping of the conventional Hanle curve onto the spatial 
domain. Thus, spins will exist primarily in the region 
where \B\ < (7T S ) _1 , i.e., where \y\ < (7T s a) _1 . The 
location of the surviving spin magnetization can be spa- 
tially shifted by adding a uniform offset to the gradient 
field. 

A measure of the resolution may be obtained from the 
width of a Gaussian fit, L r , to this curve. Fig. [3] shows 
this width as a function of the applied gradient. The 
spatial resolution L r decreases in inverse proportion to 
the gradient until L r < L s below which its rate of de- 
creases slows. This situation is reminiscent of Magnetic 
Resonance Imaging (MRI) where spatial information is 
encoded into field or frequency shifts through magnetic 
field gradients [18]. As in MRI, diffusion degrades res- 
olution, though in the current approach it can be finer 
than the spin diffusion length if one uses strong enough 
gradients, as seen in numerical simulations shown in Fig. 
[3j Furthermore, it should be noted that diffusion lengths 
for many spintronic systems are of the order of a mi- 
cron, especially at room temperatures. Field gradients 
exceeding 10 G/nm are achievable with rare-earth mi- 
cromagnetic tips [19] such as are used in Magnetic Force 
Resonance Microscopy [20-22 indicating that high spa- 
tial resolution is possible with this technique. In the case 
of imaging with a dipole tip, the gradients (and the re- 



gions of suppressed magnetization shown in Fig. [2| can- 
not be smaller than the probe-sample separation, z p: so 
for large z p this will set the resolution. 

Understanding the impact of microscopic inhomo- 
geneities on spin transport is a centrally important issue 
that calls for spatially resolved studies. Particularly im- 
portant are inhomogeneities due to the stray fields arising 
from imperfect ferromagnetic injectors used for electrical 
injection p~3j [15j [23] [24]. Spin lifetime is typically ob- 
tained from the width of Hanle curves, but such global 
measures of spin polarization are sensitive to inhomo- 
geneities, so microscopic approaches are needed to dis- 
cern the various mechanisms that can affect the shape of 
a Hanle curve. 

Non-uniformity of the interface and complex domain 
structures can cause the injected spin population to 
experience significantly different magnetic fields, and 
lose spin-information on time scales much shorter than 
spin lifetime due to dephasing. Our numerical method 
is capable of simulating the impact of magnetic fields 
with random spatial variation. Fig. [4] shows calculated 
Hanle curves in the presence of random magnetic fields. 
We consider a field B r = Af x (0, B v )x + ftf y (0,B v )y + 
A4(0, B v )z where we sample normal distributions of 
mean zero and variance B v for each field component at 
every spatial point. Fig. [4] shows results for three values 
of B v :<C (7^) _1 , = (77"s)"\ and > (7^) _1 . We see 
that the maxima can be shifted and, as in the last case, 
the curves can be significantly broadened. 

It is instructive to consider, in turn, the effects of each 
magnetic field component on the Hanle curve, as is done 
in Fig. [5| For an applied Hanle field B^ = B^y, the na ^" 
width at half-maximum £?i/2 of a Hanle curve in the pres- 
ence of an additional uniform field B n = B x x+B y y+B z z 
is given by jr s B 1/2 = y/l + ~f 2 r 2 (B x + B z ) 2 . A trans- 
verse field B x reduces S z , g (Bh = 0) (blue curve) and 
Bh ~ B x is required to significantly reduce S Zy9 further. 
A non-zero B z (red curve) on the other hand will not 
reduce S Zjg (Bh = 0) but will broaden the Hanle curve, 
since B^ ~ B z is required to significantly increase the 
spin precession cone angle. A field B y (parallel to Bh) 
shifts the peak of the Hanle curve, since S z ^ g is maxi- 
mized when the total transverse field is zero (black curve) . 
The results of Fig. [4] can be considered a superposition 
of the shifting and broadening seen in Fig. [5| Clearly 
these effects can confound the determination of spin life- 
times from Hanle widths. The detailed features of the 
stray field generated by a rough ferromagnetic injector 
will depend on several characteristics including rough- 
ness, saturation magnetization, domain size, and thick- 
ness. Regardless of these details, stray fields of the order 
kiloGauss are readily achievable a few nanometers away 
from a rough injector and these can substantially affect 
the spin polarization. In fact this has been experimen- 
tally observed in silicon devices [24] . 

In summary, we have proposed a new technique for 
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lifetime measurements, as in Hanle curves. Microscopic 
imaging of spin properties of real-world devices will be 
essential for improving their spin fidelity and functional- 
ity. 
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FIG. 4: [Color online] Global Hanle curves for a sample with 
randomly spatially varying magnetic field, B r = Af x (0, 2? v )x+ 
J\f y (0, B v )y + A/" z (0, B v )z. Such fields may be expected in 
electrical injection devices with rough injectors. The injected 
spins are assumed to be along the x-axis in this simulation. 
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FIG. 5: [Color online] Global Hanle curves for various cases of 
spatially uniform magnetic fields, B w , applied in addition to 
the swept Hanle field, Bh- These curves are valid for all values 
of D s since globally averaged spin density is insensitive to 
diffusion (see Supplementary Information), when the sample 
dimensions are much greater than L s . 



imaging spin properties of spatially inhomogeneous sam- 
ples that achieves applicability to a wide variety of ma- 
terials by relying on proven spin polarization detection 
techniques such as electrical and optical detection. The 
localized information is obtained by detecting the modifi- 
cation to the local spin density by the confined magnetic 
field of a magnetic dipole scanned over the sample. To 
this end, we have developed a method for simulating spin 
density in a medium with spatial or temporal inhomo- 
geneities of the magnetic field, lifetime, or diffusion con- 
stant, regardless of injection or detection technique. As 
in MRI, this technique has a spatial resolution inversely 
proportional to the gradient of the magnetic field down to 
length scales comparable to the spin diffusion length. Be- 
low this, resolution improves more slowly with increasing 
gradient. Our numerical analysis of spin transport em- 
phasizes the importance of experimental characterization 
of microscopic inhomogeneity, and of models that incor- 
porate these phenomena. In particular, inhomogeneities 
can confound the interpretation of spin polarization and 
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Spin dynamics in nonmagnetic materials are described 



by 



| = D s V 2 S + 7 BxS--+G, 
at t s 



(2) 



where S is the spin density, D s is the spin diffusion con- 
stant, 7 = gfiB/h is the gyromagnetic ratio, B is the 
total magnetic field experienced by the spins, r s is the 
spin relaxation time, and G represents the spin genera- 
tion term e.g., irradiation with circularly polarized light. 
S is a function of time t and spatial position r = (x,y) 
for the 2D sample that we consider. While all of the 
terms can be functions of t and r, we will consider B, G 
and r s which are independent of t, and D s will be in- 
dependent of both space and time. The formalism can 
be readily extended to include electric fields, spin-orbit 
effects, hyperfme coupling and a third spatial dimension 
for a more comprehensive treatment. Also, we are usu- 
ally interested in the steady-state solution dS/dt = and 
SV the parallel component of S. We will refer to "paral- 
lel" and "perpendicular" with respect to the injected spin 
direction. Experimentally S\\ is the most commonly mea- 
sured quantity. However, it should be noted that we can 
evaluate any component of the spin in our simulations. 

If B and r s are position-independent, the steady-state 
differential equation can be solved analytically by Fourier 
transform. In the steady state, the Fourier-transformed 
equation for the spin density S(k) takes the form 

-k 2 D s S(k) + 7 B x S(k) - S(k)/r s + G(k) = 

S(k) = [{k 2 D a + rf 1 )/ - BT^k) (3) 

where S(k) and G(k) are vectors. I is the 3x3 unit 
matrix. We introduce the 3x3 matrix 



B 



The real-space spin density is then given by the inverse 
Fourier transform of S(k). 
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An important case experimentally is given by the glob- 
ally averaged spin polarization 



Sn 



r s j A G\\{v)dh 



(4) 



for G = G\\ || and the area of integration, A, extends over 
the entire space. Note that S Zj9 oc S\\(k = 0) and thus 
from Eq. [3] is independent of D s and is given by 



S\\n = 



1 



1 



where ^ is an effective dephasing factor given by 



o2 
1 B 



1 + 72^2 



(5) 



(6) 



Spins precess in a cone whose opening half-angle is that 
between the total vector magnetic field and the injected 
spin orientation. In combination with this precession, 
the continuous injection of spins causes a distribution of 
phases, relative to ||, weighted by the spin lifetime. A 
parallel field reduces the cone opening angle resulting in 
a spin ensemble more aligned with the injected spin di- 
rection, while B± has the opposite effect. The dephasing 
factor Ob (representing the competition between parallel 
and transverse fields) describes the projection of spins 
along the injection axis and the behavior of S\\ i9 . 

The experimentally relevant situation will usually in- 
volve spatially varying fields. If B is position-dependent 
and D s ^ 0, Eq. [2] cannot be solved analytically in 
real or Fourier space for steady-state. We have found 
that solving the original time-dependent equation numer- 
ically using an Euler method works well for this problem. 
In this approach, we start from some initial condition, 
S(r, t = 0), then iterate in time using a time step At, 
according to the relation 



S(r,t + At) = S(r,t) + At x 



D s V 2 S(r,t) 



7 B(r)xS(r,f)-^^ + G(r) 



The Laplacian is evaluated numerically on a spatial grid 
of points separated by a suitable distance Ax, and is 
approximated as 

V 2 S(r, t) ~ ^ £[S(r + AxS, t) - S(r, t)] (7) 



where 5 represents the unit vectors along the spatial 
grid. We iterate until S does not change appreciably 
with time. We have verified that the analytical results 
from the analytical expressions derived above are repro- 
duced by this method. Several experimental global and 
spatially localized Hanle curves measured under varying 
conditions [3 [6] can also be qualitatively reproduced by 
this method. 
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